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Abstract 

In this article, we study the behaviour of discrete one-dimensional dynam- 
ical systems associated to functions on finite sets. We formalise the global 
orbit pattern formed by all the periodic orbits (gop) as the ordered set of pe- 
riods when the initial value thumbs the finite set in the increasing order. We 
are able to predict, using closed formulas, the number of gop for the set 
of all the functions on X. We also explore by computer experiments special 
subsets of in which interesting patterns of gop are found. 



1 Introduction 

Dynamical systems = f\{xn) associated to a function f\ : [a, 6] [a, 6] with 
(a, h) G on a real interval are very well-known. Their periodic orbits, if any, 



p 

n+l 



fx{xn) are an important feature. There is a growing interest in numerical 
analysis and industrial mathematics to study such systems, or systems in higher 
dimensions [S]. Very often, dynamical systems in several dimensions are obtained 
coupling one-dimensional ones and their properties are strongly linked. 

Simple dynamical systems often involve periodic motion. Quasi-periodic or 
chaotic motion is frequently present in more complicated dynamical systems. The 
most famous theorem in this field of research is the Sharkovskii's theorem, which 
addresses the existence of periodic orbits of continuous maps of the real line into 
itself. This theorem was once proved toward the year 1962 and published only two 
years after, A.N. Sharkovskii [TO] . 
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Mathematical results concerning periodic orbits are often obtained using sets 
of real numbers. However, most of the time, as the complex behaviour of chaotic 
dynamical systems is not explicitly tractable, mathematicians have recourse to com- 
puter simulations. The main question which arises then is: does these numerical 
computations are reliable? 

As an example O.E. Lanford III [1] reports the results of some computer exper- 
iments on the orbit structure of the discrete maps on a finite set which arise when 
an expanding map of the circle is iterated "naively" on the computer. 
Due to the discrete nature of floating points used by computers, there is a huge gap 
between these results and the theoretical results obtained when this map is consid- 
ered on a real interval. (This gap can be narrowed in some sense (i.e. avoiding the 
collapse of periodic orbits) in higher dimensions when ultra weak coupling is used 

mm 

Nowadays the claim is to understand precisely which periodic orbit can be ob- 
served numerically in such systems. In a first attempt we study in this paper the 
orbits generated by the iterations of a one-dimensional system on a finite set X^v 
with a cardinal N. The final goal of a good understanding of the actual behaviour of 
dynamical systems acting on fioating numbers (i.e. the numbers used by computers) 
will be only reached when this first step will be achieved. 

On finite set, only periodic orbits can exist. For a given function we can compute 
all the orbits, all together they form a global orbit pattern. We formalise such a 
gop as the ordered set of periods when the initial value thumbs the finite set in the 
increasing order. We are able to predict, using closed formulas, the number of gop 
for the set jF/y of all the functions on X. We also explore by computer experiments 
special subsets of J-'n in which interesting patterns of gop are found. 

This article is organized as follows: in section [2] we display some examples already 
known of such computational discrepancies for the logistic and circle maps. In 
section [3] we introduce a new mathematical tool: the global orbit pattern, in order 
to describe more precisely the behaviour of dynamical systems on finite sets. In 
section m we give some closed formulas giving the cardinal of classes of gop of J-'n- 
In section |5] we study the case of sets of functions with local properties versus their 
gop, in order to show the significance of these new tools. 
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2 Computational discrepancies 

2.1 Approximated logistic map 

As an example of such collapsing effects, O.E. Lanford III, presents the results of a 
sampling study in double precision of a discretization of the logistic map 

x^+i = 1-2x1 (1) 

which has excellent ergodic properties. The precise discretization studied is obtained 
by first exploiting evenness to fold the interval [—1, 0] to [0, 1], i.e replacing ([T]) by 

x„+i = |l-2a;2| (2) 

On [0, 1] it is not difficult to see that the folded map has the same periods as the 
original one. The working interval is then shifted from [0, 1] to [1, 2] by translation, 
in order to avoid perturbation in numerical experiment, by the special value 0. Then 
the translated folded map is programmed in straightforward way. 
Out of 1,000 randomly chosen initial points, 

• 890, i.e., the overwhelming majority, converged to the fixed point correspond- 
ing to the fixed point-1 in the original representation ([T]), 

• 108 converged to a cycle of period 3,490,273, 

• the remaining 2 converged to a cycle of period 1,107,319. 

Thus, in this case at least, the very long-term behaviour of numerical orbits is, 
for a substantial fraction of initial points, in fiagrant disagreement with the true 
behaviour of typical orbits of the original smooth logistic map. 

2.2 Circle approximated maps 

In the same paper, O.E. Lanford III, studies very carefully the numerical approxi- 
mation of the map 

Xn+i = 2x„ + 0.5x„(l — Xn) (mod 1) < X < 1 (3) 

It is perhaps better to think this map (see Figured]) as acting on the unit interval 
with endpoints identified, i.e., on the circle. Note that f'{x) > 1.5 everywhere, so / 
is strictly expanding in a particularly clean and simple sense. As a consequence of 
expansivity, this mapping has about the best imaginable ergodic properties : 
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Figure 1: Graph of the map f{x) = 2x + 0.5x(l — x) (mod 1) 

• it admits a unique invariant measure /i equivalent to Lebesgue measure, 

• the abstract dynamical system (/, fi) is ergodic and in fact isomorphic to a 
Bernoulli shift, 

• a central limit theorem holds. 

One consequence of ergodicity of / relative to fi is that, for Lebesgue almost 
initial points x in the unit interval, the corresponding orbit /*•"''* (x) is asymptotically 
distributed over the unit interval according to fi. 

In numerical experiments performed by the author, the computer working with 
fixed finite precision is able to represent finitely many points in the interval in 
question. It is probably good, for purposes of orientation, to think of the case where 
the representable points are uniformly spaced in the interval. The true smooth map 
is then approximated by a discretized map, sending the finite set of representable 
points in the interval to itself. 

Describing the discretized mapping exactly is usually complicated, but it is 
roughly the mapping obtained by applying the exact smooth mapping to each of 
the discrete representable points and "rounding" the result to the nearest repre- 
sentable point. 

O.E. Lanford III has done two kinds of experiments, in the first one he uses 
double precision floating points, in the second uniformly spaced points in the interval 
with several order of discretization (ranging from 2^^ to 2^^) are involved. In each 
experiment the questions addressed are: 
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Double precision (sampling) 
7 cycles found 



period 


basin size 


relative size 


27,627,856 


517 


51.7 % 


88,201,822 


320 


32.0 % 


4,206,988 


147 


14.7 % 


4,837,566 


17 


1.7% 


802,279 


8 


0.8 % 


6,945,337 


6 


0.6 % 


2,808;977 


1 


0.1 % 



Table 1: Coexisting periodic orbits of mapping (jl]) for double precision numbers. 

• how many periodic cycles are there and what are their periods ? 

• how large are their respective basins of attraction, i.e., for each periodic cy- 
cle, how many initial points give orbits with eventually land on the cycle in 
question? 

On one hand, for the experiments using ordinary (IEEE- 754) double precision - 
so that the working interval contains of the order of 10^^ representable points - 1,000 
initial points at random are used in order to sample the orbit structure, determining 
the cycles to which they converge. 

The computations were accomplished by shifting the working interval from the 
initial interval [0, 1] to [1, 2] by translation, for the same previous reason as for the 
logistic map, i.e., the map actually iterated was 

Xn+i = 2x„ + 0.5x„(a;„ — 1)(2 — a;„) (mod 1) 1 < x < 2 (4) 
Some results are displayed in Table [H 

On another hand, for relatively coarse discretizations the orbit structure is de- 
termined completely, i.e., all the periodic cycles and the exact sizes of their basins 
of attraction are found. Some representative results are given in Tables [2] to [51 In 
theses tables, N denotes the order of the discretization, i.e., the representable points 
are the numbers, with < j < A^. 

It seems that the computation of numerical approximations of the periodic orbits 
leads to unpredictable results. 

Many more examples could be given, but those given may serve to illustrate the 
intriguing character of the results: the outcomes prove to be extremely sensitive to 
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441 4S9 


64 85 % 


1,746 


2,946,734 


35.13 % 


13 


205 


24.10-^ % 


6 


132 


16.10-^ % 


30 


96 


11.10-4 % 


4 


8 


< 1.10-6 % 


1 


1 


< 1.10-6 % 



Table 2: Coexisting periodic orbits of mapping (jll) for the discretization N = 2' 

N = 2'^^ = 16, 777, 216 
2 cycles 

period basin size relative size 
5,300 16,777,214 100 % 

1 2 < 1.10-6 % 



Table 3: Coexisting periodic orbits of mapping (j4]) for the discretization = 2^'*. 

the details of the experiment, but the results all have a similar flavour: a relatively 
small number of cycles attracts near all orbits, and the lengths of these significant 
cycles are much larger than one but much smaller than the number of representable 
points. O. E. Lanford III wrote that there are here regularities which ought to be 
understood. 

In pLj, P. Diamond and A. Pokrovskii, suggest that statistical properties of the 
phenomenon of computational collapse of discretized chaotic mapping can be mod- 
eled by random mappings with an absorbing centre. The model gives results which 
are very much in line with computational experiments and there appears to be a 
type of universality summarised by an Arcsine law. The effects are discussed with 
special reference to the family of mappings 

Xn+l = I - \l - 2Xnf 0<X <l l<i<2 (5) 

Computer experiments show close agreement with prediction of the model. 
However these results are of statistical nature, they do not give accurate infor- 
mation on the exact nature of the orbits (e.g. length of the shortest one, of the 
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N = 


224 _ 1 = 


,777,215 




10 cycles 




period 


basin size 


relative size 


3,081 


7,502,907 


44.72 % 


699 


3,047,369 


18.16 % 


3,469 


2,905,844 


17.32 % 


1,012 


2,774,926 


16.54 % 


563 


290,733 


11.73 % 


2,159 


221,294 


1.32 % 


138 


21,610 


0.13 % 


421 


12,477 


0.07 % 


9 


54 


< 1.10-3 % 


1 


1 


< 1.10-7 % 



Table 4: Coexisting periodic orbits of mapping (jlj) for the discretization = 224 — 1. 



N 


= 225 = 33, 


554,432 




8 cycles 


period 


basin size 


relative size 


4,094 


32,114,650 


95.71 % 


621 


918,519 


2.74 % 


283 


516,985 


1.54 % 


126 


2,937 


< 0.01 % 


6 


887 


< 0.01 % 


55 


433 


< 0.01 % % 


4 


20 


< 1.10-6 % 


1 


1 


< 1.10-*^ % 



Table 5: Coexisting periodic orbits of mapping (jlj) for the discretization N = 2' 
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greater one, size of their basin of attraction). It is why we consider the problem of 
computational discrepancies in an original way in the next section. 

3 Global orbit pattern 

We introduce in this section a new mathematical tool: the global orbit pattern of 
a function, in order to describe more precisely which kind of behaviour occurs in 
discretized dynamical systems on finite sets. 

3.1 Definitions 

Let / be a map from X onto X. 

We denote by Vx e X, f{x) = x and Vp > 1, Vx e X, fP{x) = / o / o^. . . o / (x). 

p times 

Given xq E X, wc define the sequence {xi} of the associated dynamical system by 

Xi+i = f{xi) for i > 0. Thus x., = p{xo). 

The orbit of xq under / is the set of points O{xo, /) = {/^(xo), i > 0}. 
The starting point Xo for the orbit is called the initial value of the orbit. 
A point X is a fixed point of the map / if f{x) — x. 

A point a; is a periodic point with period p if /^(x) = x and /'^(x) ^ x for all k such 

that < k < p, p is called the order of x. 

If X is periodic of order p, then the orbit of x under / is the finite set 

{x, f{x),p{x),..., p~^{x)}. We will call this set the periodic orbit of order p or a 

p-cycle. 

A fixed point is then a 1-cycle. 

The point x is an eventually periodic point of / with order p if there exists K > Q 
such that \/k>K f^+P{x) = f^{x). 

\/x e X, we denote a;(x, /) the order of x under / or simply a;(x) when the map / 
involved is obvious. 

A subset T of X is invariant under / if f~^{T) = T. That is equivalent to say that 
T is invariant under / if and only if /(T) C T and f-\T) C T. 

3.2 Map on finite set 

Along this paper, N is a, non-zero integer and stands for the cardinal of any finite 
set A. In this article, we consider X as an ordered finite set with elements. We 
denote it Xj^ and Xj^ is isomorphic to the interval in N [0, A^ — 1]. Then ^Xj^ — N. 
Let / be a map from X^ into X^. We denote by J^n the set of the maps from 
Xn into Xn. Clearly, J^n is a finite set and Ji^jv = A^^ elements. For all x e X^, 
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0{x, f) is necessarily a finite set with at most elements. Indeed, let us consider 
the sequence {x, f{x),p{x),..., f'^~^{x), f^{x)} of the first + 1 iterated points. 
Thanks to the Dirichlet's box principle, two elements are equals because Xn has 
exactly N different values. Thus, every initial value of leads ultimately to a 
repeating cycle. More precisely, if a; is a fixed point 0{x, /) is the unique element 
X and if a; is a periodic point with order p, 0{x, /) has exactly p elements. In 
this case, the orbit of x under / is the set 0{x, f) = {x, /(x), P{x), . . . , /^~^(a;)}. 
If X is an eventually periodic point with order p, there exists > such that 
> K f^^^{x) = f^{x). In this case, the orbit of x under / is the set 0{x, /) = 
{x, /(x), . . . , /^(x), f^+\x), . . . , f^^^-\x)}. 



3.3 Equivalence classes 
3.3.1 Components 

Let / G Tn- We consider on X^ the relation ~ defined by : Vx, x' G X^, x ~ x' -v^ 
3k en such that /^(x) G 0(x', /). The relation ~ is an equivalence relation on Xjy. 
Sn / ~ is the collection of the equivalence classes that we will call components of X^ 
under / which constitute a partition of Xjy. The number of components are given 
in |3j. Asymptotic properties of the number of cycles and components are studied in 
[8]. For each component, we take as class representative element the least element of 
the component. The components will be written T/v(xo, /),..., T/v(xp^ ^, /) where 
Xi is the least element of T^i^Xi, /). 

By analogy with real dynamical systems, we can define attractive and repulsive 
components in discretized dynamical systems as follows. 

Definition 1 A component is repulsive when it is a cycle. Otherwise, the 
component is attractive. 

Remark In other words, a component is attractive when the component con- 
tains at least an eventually periodic element. The corresponding cycle is strictly 
contained in an attractive component. 

Examples are given in the tables [6], [7] and [HI 

For instance, in the table[7l the function / has {8} as orbit and {1, 7, 8} as com- 
ponent which is attractive because {1} and {7} are eventually periodic elements. 
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Function orbit/component/nature 
0^4 
1^1 
2^5 
3^4 
4^0 
5^9 
6^8 
7^5 
8^6 
9^7 



{0,4} 


{0,4,3} 


attractive 


{1} 


{1} 


repulsive 


{5,9,7} 


{2,5,9,7} 


attractive 


{6,8} 


{6,8} 


repulsive 



Table 6: Orbits and components of a function belonging to Tio with gop [2, 1, 3, 2]io. 



Function orbit/component/nature 
0^6 
1 ^ 8 



2^4 
3^9 
4^5 
5^2 
6^0 

7 ^ 1 



{0,6} 


{0,6} 


repulsive 


{8} 


{1,7,8} 


attractive 


{2,4,5} 


{2,4,5} 


repulsive 


{3,9} 


{3,9} 


repulsive 



Tabic 7: Orbits and components of another function belonging to J^io with gop 
[2,l,3,2]io. 

3.3.2 Order of elements 

Here are some remarks on the order of elements of components. 

Remark The order of every element of a component is the length of its inner 
cycle. 
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{2,9} 


{0,2,9} 


attractive 


{6} 


{1,6} 


attractive 


{3,8,4} 


{3,8,4} 


repulsive 


{5,7} 


{5,7} 


repulsive 



Function orbit / component / nature 

0^9 

1 ^ 6 
2^9 
3^8 
4^3 
5^7 
6^6 
7^5 
8^4 
9^2 



Table 8: Orbits and components of a third function belonging to Tiq with gop 
[2,l,3,2]io. 

Definition 2 For all x G X^r, there exists i G |IO,p/^7v] such that x belongs to 
the component T^i^Xi, /). Then u^x, f) is equal to the order u^Xi, /). 

Remark For all i G |0,]5/,Ar], Tj^{xi, f) is an invariant subset of Xj^ under /. 

In the example given in the table [6], the order of the element {5} is 3, the order 
of the element {1} is 1, the order of the element {8} is 2. The elements {2} and {5} 
have the same order. 



3.4 Definition of global orbit pattern 

For each / G J^n, we can determine the components of Xn under /. For each com- 
ponent, we determine the order of any element. Thus, for each / G J-'n, we have a 
set of orders that we will denote Q{f,N). Be given /, there exist pf^N components 
and pf^N representative elements such that Xq < Xi < . . . < Xp^^. 

For each / G J-'n, the sequence [ti;(a;o), . . . ,u{xpj_. f]jr^ with xq < xi < 

. . . < Xpj will design the global orbit pattern of / G J^n- 

We will write gop{f) = [uj{xo),uj{xi), . . . ,u{xp^,^); f]jr^. 

When the reference to / G J^n is obvious, we will write shortly 
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gopif) = [uj{xo),uj{xi), . . . ,uj{xp^^^)]n or gop{f) = [uj{xo),uj{xi), . . . ,uj{xp^.^)] . 

For example, the same gop associated to the functions given in the tables O [7] 
and [8] is [2, l,3,2]io. 

Another example is given in the table O In that example, we have co'(O) = 2, 
u;(3) = 1, u;(4) = 4. 



Function orbit/component/nature 

^ 1 

1-^0 {0,1} {0,1,2} attractive 

2^0 

3-^3 {3} {3} repulsive 

5^6 {4,5,6,7} {4,5,6,7} repulsive 

6^7 

7^4 



Table 9: Orbits and components of a function belonging to jFg with gop [2, 1,4]8. 

Definition 3 The set of all the global orbit patterns of jFjy is called Q{Tn)- 

For example, for = 4, the set ^(^^4) is {[1]; [1, 1]; [1, 1, 1]; [1, 2]; [1, 1, 1, 1]; 
[1, 1, 2]; [1, 2, 1]; [1, 3]; [2]; [2, 1]; [2, 1, 1]; [2, 2]; [3]; [3, 1]; [4]}. 

3.5 Class of gop 

We give the following definitions : 

Definition 4 Let he A= [wi, . . . , uj^n a gop. Then the class of A, written A, 
is the set of all the functions / G J^n such that the global orbit pattern associated 
to / is A. 



For example, for N = 10, the class of the gop [2, 1, 3, 2]-^q contains the following 
few of many functions defined in the tables [6l [7] and [HI The periodic orbit which are 
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encountered have the same length nevertheless there are different. 

Definition 5 Let he A — [uji, . . . ,u;p\n a gop. 

p 

Then the modulus of A is \A\ — X] 

fc=i 

Remark | [a;i, . . . , a;^] jv| < N. 



3.6 Threshold functions 

3.6.1 Ordering the discrete maps 

Theorem 1 The sets J^n and |1, N^j are isomorphic. 

Proof. We define the function from J-'n to |1, A^^] by : for each / G J-'n, 4>{f) is 

N-l 

the integer n such that n = ^ f{k)N^~^~'' + 1. 

k=0 

Then is well defined because n G |1, N^]. 

Let n be a given integer between 1 and A^^. We convert n — 1 in base N : 
there exists a unique iV-tuple (a„_i,o; cin-i.i', ■ ■ ■', cin-i,N-i) G |0, iV — 1]^ such that 

N-l 

n — 1 an-i^N-i-iN^~^~^ . We can thus define the map with : Vi e X^v, 

i=0 

= an-i,N-i-i- Then is one to one. 
Remark This implies J-'n is totally ordered. 

Definition 6 Let / e J^^v- Then n = ^ f{k)N^-'^-^ + 1 is called the rank 

fe=0 

of/. 



3.6.2 Threshold functions 

Be given a global orbit pattern A, we are exploring the class A. 

Theorem 2 For every A e Q{!Fn), the class A has a unique function with 
minimal rank. 

Definition 7 For every class A e Q{Tn)-, the function defined by the previous 
theorem will be called the threshold function for the class A and will be denoted by 
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Tr{A) or Tr{A). 

To prove the theorem, we need the following definition : 

Definition 8 Let / G J-'n be a function. Let p be a non zero integer smaller 
than N. Let he Xi, . . . ,Xp p consecutive elements of Xj^. Then xi, . . . , is a canon- 
ical p-cycle in relation to / if Vj G — 1] , f{xj) = Xj+i and f{xp) = xi. 

Proof. Let [ui, . . . ,ujp] be a global orbit pattern of Q{J^n)- We construct a specific 
function / belonging to the class [ui, . . . , Up] and we prove that the function so ob- 
tained is the smallest with respect to the order on ^jv- With the first Ui elements of 
|0, — 1], that is the set of integers [0, cji — 1], we construct the canonical Wi-cycle: 
if ui = 1, we define /(O) = 0, else /(O) = 1, /(I) = 2, . . ., /(tui - 2) = c^i - 1, 
/(^i - 1) = 0. 

Then Vj G [cji - 1, cji + - s - 1], we define f{j) = 0. 

Then with the next UJ2 integers [ui + N — s,u!i + N — s + 002 — l] we construct the 
canonical u;2-cycle. We keep going on constructing for all j G |3,p] the canonical 
Uj-cjde. 

In consequence, we have found a function / belonging to the class [ui, . . . ,u!p]. 
Assume there exists a function g G J-'n belonging to the class of / such that g < f ■ 
Let J = G [0, A^ — 1] such that f{i) 7^ 0}. As g < f, there exists io E I such that 
fi'(^o) < f{io)- There exists also jo such that io G uj^. If f{io) = io, then = 1, 
g{io) < io and then g{io) ^ ujj^. Then the global orbit pattern of g doesn't contain 
anymore 1 as cycle. The global orbit pattern of g is different from the global orbit 
pattern of /. If f{io) = io + l, then g{io) < io- Either g{io) = io and then the global 
orbit pattern of g is changed, or g{io) < io and we are in the same situation as pre- 
viously. Thus, in any case, the smallest function belonging to the class [ui, . . . ,Up] 
is the one constructed in the first part of the proof. 

The proof of the theorem gives an algorithm of construction of the threshold 
function associated to a given gop. 

The threshold function associated to the gop [2, l,3,2]io is explained in the table 
M Its rank is n = 1, 000, 467, 598. 

Theorem 3 There are exactly 2^ — 1 different global orbit patterns in ^jy- 

That is tl^(J^jv) =2^-1. 



For example, for A^ = 4, tl^(J^4) = 2^ - 1 = 15. 



126 



Rene Lozi, Clarisse Fiol 



remaining 
images 



First step Second step Third step Fourth step Fifth step 

(Jonstruction Construction Construction (Jonstruction Filling the 

of the first of the last of the of the 

canonical canonical canonical canonical 

2-cycle 2-cycle 3-cycle 1-cycle with 

0^1 0^1 0^1 0^1 0^1 

1^0 1^0 1^0 1^0 1^0 

2^ 2^ 2^ 2^ 2^0 

3^ 3^ 3^ 3^ 3^0 

4^ 4^ 4^ 4^4 4^4 

5^ 5^ 5^6 5^6 5^6 

6^ 6^ 6^7 6^7 6^7 

7^ 7^ 7^5 7^5 7-^5 

8^ 8^9 8^9 8^9 8^9 



Tabic 10: Algorithm for the threshold function construction for the gop [2, 1, 3, 2]io. 



Proof. Let p be an integer between 1 and N. Consider the set L{p, N) of p-tuples 

(«!, . . . , Op) e (N*y such that ai + . . . + ap < N. 

We write L{N) — {L{p, N),p — 1 . . . N}. L{N) and Q{!Fn) have the same elements. 
Then 

p=N p=N , . 

p=l p=l \ / 

3.6.3 Ordering the global orbit patterns 

We define an order relation on Q{Tn)- 

Proposition 1 Let A and B be two global orbit patterns of Q (J^n) ■ 
We define the relation -< on the set Q{J^n) by 

A^B iff Tr{A) < Tr{B) 
Then the set {Q{Tn)-, -<) is totally ordered. 
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Proof. As the order -< refers to the natural order of N, the proof is obvious. 

Let r > 1, p > 1 be two integers. Let [ui, . . . ,Ljp] and . . . ,ul.] be two 
global orbit patterns of Q{Tn)- For example, if p < r, in order to compare them, 
we admit that we can fill [cui, . . . ,ujp\ with r — p zeros and write [ui, . . . ,u!p] = 
[ui, . . . ,t^p,0, . . . ,0]. 

Proposition 2 Let r > 1, p > 1 be two integers such that p < r. Let 
A = [uji, . . . ,ujp] and B = [uj[, . . . , ujI] be two global orbit patterns. 

• If r = p = 1 and uji < uj[ then A -< B. 

• If r > 2 then 

* If cji < uj[ then A -< B. 

* If cui = uj[ then there exists K G [2;r] such that ujk 7^ ^jJ'k ci-nd 
Vi < K, LUi = uo[. 

• If \A\ < \B\, then A^B. 

• If l^l = \B\, then if cok < ^'k then A -< B. 

For example, for = 4, the global orbit patterns are in increasing order : 
[1] ^ [1, 1] ^ [1, 1, 1] ^ [1, 2] -< [1, 1, 1, 1] ^[1,1,2]^ [1, 2, 1] -< [1, 3] -< [2] -< [2, 1] -< 
[2, 1,1] ^[2, 2] ^[3] ^[3,1] ^[4]. 

3.6.4 Algorithm for ordering the global orbit patterns : a pseudo-decimal 
order 

The table [H] gives a method for ordering the gop : indeed, we consider each gop as 
if each one represents a decimal number: we begin to order them in considering the 
first order uJi. Considering two gops A= [tui, cjp] and A' = [uj[, ...,uj'^], if oji < uj[, 
then A -< A'. For example, [1, 3] -< [2, 1]. If uji = uj[ and \A\ — ui < \A'\ — uj[, then 
A -< A'. For example to compare the gop [1, 2] and the gop [1, 1, 1, 1], we say that 
the first order Ui stands for the unit digit - which is cui = 1 here, then the decimal 
digits are respectively 0.2 and 0.111. We calculate for each of them the modulus-a;i: 
we find |[1,2]| - 1 = 2 and |[1,1,1, 1]| -1=3, thus [1,2] -< [1,1,1,1]. Finally, if 
uji = uj[ and \A\ — ui = \A'\ — uj'i, then also we use the order of the decimal part. 
For example, [1, 1, 1] -< [1, 2] because 1.11 < 1.2. Applying this process, we have the 
sequence of the ordered gop for A^ = 4 given in the previous paragraph. 

For example, for A^ = 4, we construct a tree (Figure [2]) : each vertex is an ordered 
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[1,1,1,1] 



Figure 2: Tree for the construction of the order on Q{J^i) 
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rinn 


iVlUQUlUb 


iVi U U. U 1 U o- Uy 1 




iViUU-UlUo 


iVi U U. U. 1 U o- U/ 1 


[-■-J 


1 
1 


n 

U 




9 


n 
u 


fl 11 




1 
± 


\2 11 




1 

X 




3 


2 


[2, 1, 1] 


4 


2 


[1,2] 


3 


2 


[2,2] 


4 


2 


[1,1,1,1] 


4 


3 








[1,1,2] 


4 


3 


[3] 


3 





[1,2,1] 


4 


3 


[3, 1] 


4 


1 


[1,3] 


4 


3 


[4] 


4 






Table 11: Ordered gop for = 4 with modulus and modulus-o;! 



orbit, the modulus of the gop is written on the last edge. 

However, the sequence of ordered gop differs from the natural downward lecture of 
the tree and has to be done following the algorithm. 



4 Numbering the discrete maps 

In this section we emphasize some closed formulas giving the cardinal of classes of 
gop. Recalling first the already known formula for the class [ 1, . , IJ at for which 

k times 

we give a detailed proof, we consider the case were the class possesses exactly one 
fc-cycle, the case with only two cycles belonging to the class and finally the main 
general formula of any cycles with any length. We give rigorous proof of all. The 
general formula is very interesting in the sense that even using computer network 
it is impossible to check every function of jF/v when is larger than 100 or may be 
50. 



4.1 Discrete maps with 1-cycle only 

The theorem 4 gives the number of discrete maps of J-'n which have only fixed points 
and no cycles of length greater than one. This formula is explicit in ^ and [9j. A 
complete proof is given here in detail. 

Theorem 4 Let k be an integer between 1 and A^. The number of functions 
whose global orbit pattern is [ 1, . ^. , Ij jy (i.e. belonging to the class [ 1, . ^. , Ij jy) is 

k times k times 
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N-1 
N-k 



N-k 



That is ^\1^_^N 

k times 



N-1 
N-k 



N 



N-k 



Proof. Let A; be a non-zero integer. Let / be a function of Tjy. There are 



N 
k 



possibihties to choose k fixed points. There remain N-k points. Let p be an integer 
between 1 and N-k. We assume that p points are directly connected to the k fixed 
points. For each of them, there are k manners to choose one fixed point. There are 
k^ ways to connect directly p points to k fixed points. There remains N — k—p points 
that we must connect to the p points. There are ji[l, . . . , Ijjv-fe functions. Finally, 



N 



the number of functions with k fixed points is J k^ ti[ l, . ^. , ij jy-fc- We now 

^ ^ n times 



p times 
N-k 



prove recursively on N for every < k < N that |j| l, , I Jat 

k times 

We have ji[l]^ = 1. The formula is true. 

'N-1 
N-k 



p times 

N-1 
N-k 



N 



N-k 



We suppose that yk < N W[l,...,l\ 



N 



N-k 



k times 

Let X be a set with + 1 elements. Wc look for the functions of J^n+i which have 
k fixed points. Thanks to the previous reasoning, we have 



k times 



^\1^_-A}n+i = 

k times 

We use the recursion assumption. 



N+1 
k 



N+1 
k 

N+1 
k 



N+l-k 



p=l 

N+l-k 



[ 1, ■ - , Ij Af+l-fc- 
p times 



J2 F tt[ir^iv-(fe-i). 

p times 



[ 1; ) IJ aT+I 
k times 



N+l-k 



E k^[1_';){N-k+i)--''-^^-^. 




k 



p 

N-k 



k 
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k times 



N+1 



N-k 



k times 

N+1 



k times 



N 

N-k 



4.2 Discrete maps with /c-cycle 

We look now for the number of functions with exactly one /c-cycle. 

Theorem 5 Let k be an integer between 1 and N. The number of functions 
whose global orbit pattern is [k]N is tl[ li ■ - ; IJ a^ x (A; — 1)!. 

k times 



i.e. mN = ^^h^^NX{k~ 1)1 



k times 



Proof. There are (^) ways of choosing k elements among N. Then, there are {k — l)\ 
choices for the image of those k elements in order to constitute a /c-cycle by /. We 
must now count the number of ways of connecting directly or not the remaining 
N — k elements to the /c-cycle. We established already this number which is equal 

^ _ N-k 

to J2 tl[ir^iv-fc. Finally, we have mN ^ (k - E ttU^-^iv-fe. 

^ p times p times 

That is, mN = ^\1^_^N X (A; - 1)!. q.e.d 

k times 

4.3 Discrete maps with only two cycles 

We give the number of functions with only two cycles. 

Theorem 6 Let N > 2. Let p and q be two non-zero integers such that 
p + q< N. Then, 

ttb, qU = Ul^_^\N^^ = {N-{p+g)y. q ■ 
p+q times 
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Proof. We consider a function / which belongs to the class [ 1, . . . , 1 Jat. We search 

p+q times 

the number of functions constructed from / whose gop is [p, q]^- From the p fixed 
points of /, we construct a p-cycle. Thus, there arc {^^1^^^) ways to choose p — 1 
integers among the p + q — 1 fixed points. Counting the first given fixed point of /, 
we have p points which allow to construct (p — 1)! functions with a p-cycle. Then 
there remain q points which give {q — 1)! different functions with a g'-cycle. Finally, 
the number of functions whose gop is [p, q]^ is: — l)\{q — 1)1 that is the 

formula (ttlll. 

Q 

Remark We notice that for all k non-zero integer such that k < N — 1, 



4.4 Generalization : discrete maps with cycles of any length 

We introduce now the main theorem of the section which gives the number of gop 
of discrete maps thanks to a closed formula. 



Given a global orbit pattern a, the next theorem gives a formula which gives the 
number of functions which belong to a. 



Theorem 7 Let p > 2 be an integer. Let [cui, . . . ,u!p]N be a gop of Q{J-n)- 

Then, 

{oji+...+u}p—iy. 



PiN 



N 



U)pX{ojp-l+ljJp)x...x{u>2 + ---+Ulp) 



u}i+...+u!p times 

(TV— 1)! Ar-'^~('^i+---+'^p' 
{N-{m+...+u,p)y. n(E^i) 

k=2 j=k 



Proof. We consider a function / which belongs to [1, , lljv. We search 



uj\+...+LOp times 

the number of functions constructed from / whose gop is . . . ,oJp\N- From the 
uji fixed points of /, we construct a a;i-cycle. Thus, there are ways to 

choose cji — 1 integers among the uoi -\- . . . + ujp — 1 fixed points. Counting the first 
given fixed point of /, we have uii points which allow to construct {uji — 1)! functions 
with a a;i-cycle. Then, the first fixed point of / which has not be chosen for the 
a;i-cycle, will belong to the a;2-cycle. Thus, there are ('^^^^^^-T"^) ways to choose 
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a;2 — 1 integers among the 002 + ■ ■ ■ + (jJp — l fixed points. So we have 002 points which 
allow to construct {002 — 1)' functions with a a;2-cycle. We keep going on that way 
until there remain ujp fixed points which allow to construct {cUp — 1)! functions with 
a a;p-cycle. Finally, we have constructed : 

r^^;t?-^)(a;i-l)!r+^,+?-^)(a;2-l)! x . . . x (--;+?-^) - 1)!^- 1)! func- 
tions. We simplify and obtain the formula. 

Corollary 1 Let p be a non-zero integer. Let [c^-'i. . . . , ujp]N be a gop of Q{J-n). 
We suppose that there exists j such that cuj >2. Let h be an integer between 1 and 
ujj — 1. Then 



tl[a;i, . . .,ujj, . . .,uJp]j^ = ^[uJi, . . .,ujj - h,h,ujj+i, . . . ,0;^]^ x {h + ujj+i + . . .+ujp). 



Proof. - h,h,u}j+i, . . . ,Up]j^ X {h+u}j+i + . . .+u}p) = (ifl, ,1 



N 



X 



u)i+...+u}p times 

(uii+...+uip-l)\(h+ujj+i+...+ujp) 



(jjp{ujp-i+ujp)...{ojj+i+...+(jjp){h+ujj+i+...+ujp){wj+u>j+i+...+ujp)x...x ' 

We simplify and we exactly obtain 



...,Uj-h,h, Uj+i, Up]j^ X {h + Uj+i + . . .+Up) = ^[ui, ...,u}j,..., Up]j^. 
Examples : 



ii[2,l,3,2]io = 302,400. 

tt[5,2, 10,8, 15,2, 3]5o = 29, 775, 702, 147, 667, 389, 218, 762, 343, 520, 975, 006, 348, 329, 
578, 044, 480, 000, 000, 000, 000, 000. 



tt[5, 2, 10, 8, 15, 2, 3)50 ^ 2.98 x 10^^ among the 8.88 x 10^^ functions of 



5 Functions with local properties 
5.1 Functions with locally bounded range 

Since several centuries, continuity and differentiability play a dramatic role in math- 
ematical analysis. However these concepts are not transposable to the functions on 
finite sets. Then in this section, in order to obtain more precise results on the orbit 
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of dynamical systems on finite sets, we introduce subsets of jF^r, whose functions 
have local properties such as locally bounded range. In these subsets, the gop are 
found to be fully efficient in order to describe very precisely the dynamics of the 
orbits. We first consider the very simple subset Ci^n of functions for which the 
difference between f{p) and f{p+ 1) is drastically bounded. In subsection 15.21 we 
consider more sophisticated subsets. 

We consider the set : 

£i,7v = {feJ^N such that Vp, < p < iV - 2, \f{p) - /(p + 1) | < 1}. 



5.1.1 Orbits of £i at 

Theorem 8 If / G Ci^n then / has only periodic orbits of order 1 or 2. 

Proof. We suppose that / G Ci^n has a 3-cycle. We denote {a; f{a); f'^{a)) taking 
a the smallest value of the 3-cycle. If a < /(a) < P{a) then there exist two 
non-zero integers e and e' such that /(a) = a + e and P{a) = f{a) + e'. Thus, 
/2(a) - e' < f{a) < p{a) + e'. That is /(a) <a< f{a) + 2e'. And finally we have 
the relation a + e < a which is impossible. 

If a < P{a) < f{a) then there exist two non-zero integers e and e' such that 
/2(a) = a + e and /(a) = p{a) + e'. Thus, /(a) - e < f{a) < f{a) + e. That is 
/(a) — e < a < f{a) + e. But /(a) — e = a + e'. And finally we have the relation 
a + e' < a which is impossible. 

We can prove in the same way that the function / can't have either 3-cycle or greater 
order cycle than 3. 

5.1.2 Numerical results and conjectures 

We have done numerical studies of the Q{Ci^n) for = 1 to 16, using the brute 
force of a desktop computer (i.e. checking every function belonging to these sets). 

The tables [121 [131 [HI [13 [HI and [T7| show the sequences for Ci^i to >Ci,i6- 

In theses tables we display in the first column all the gop of Q{Ci^n) for every 
value of A^. For a given A^, there are two columns; the left one displays the cardinal 
of every existing class of gop (- stands for non existing gop). Instead the second 
shows more regularity, displaying on the row of the gop [2,2, ...,2] the sum of the 




k times 
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cardinals of all the classes of the gop of the form [2,2, , 1 , . . . , 2] which exist. 



th 



I 

V 

fc+1 orders 



Then we are able to formulate some statements which have not yet been proved. 

1 if A; = 1 



Statement 1 Jj c^ n = S ifA; = 2 

1^^;::^^^ ' I (^)(^ + l)x3^ for 3<k<^ 



Remark We call Uk — 1, , ^jciN- ^^''^ k > 2, then Uk is the sequence 



N—k+l times 

A120926 On-line Encyclopedia of integer Sequences : it is the number of sequences 
where is isolated in ternary words of length written with {0, 1, 2}. 



Statement 2 ji 


11. 


1 








k times 




Statement 3 ft 


12- 


2 


'2J 



£-l.N 



1,1, ,1 



for k < ^. 



k+l times 



,N+2 



N 
2 ■ 



k times 



times 



Statement 4 



ttL2,2, ,2Ja„ = tt[2,2, ,2,lj£,_ for 2A; < TV < 3A; - 1. 



k times 

Statement 5 



k times 



k+1 



2,2, ,2\c, , = E tt[2,2, ,J.^,...,2\c,^, for 2A; + 1 < TV. 

i=i 



k times 



i 

fc+l orders 



These statements show that first the set >Ci^jv is an interesting set to be consid- 
ered for dynamical systems and secondly the gop are fruitful in this study. However 
the set 



C2,N ^{fe^N such that Vp, < p < TV - 2, - /(p + 1) | < 2} 

is too much large to give comparable results. Then we introduce more sophisti- 
cated sets we call sets with locally bounded range which more or less correspond to 
an analogue of the discrete convolution product of the local variation of / with a 
compact support function a^. 
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g.o.p. 


N=l N=l 


N=2 


N=2 


N=3 


N=3 


N=4 


N=4 


Total number 


1 




4 




17 




68 


[1] 


1 + 


2 


+ 


7 


+ 


26 


+ 


[1,1] 


+ 


1 


+ 


4 


+ 


14 


+ 


[1,1,1] 


+ 




+ 


1 


+ 


4 


+ 


[1,1,1,1] 


+ 




+ 




+ 


1 


+ 


[2] 


+ 


1 


1 


4 


4 


18 


18 


[2,1] 


+ 




+ 


1 


1 


3 


4 


[1,2] 


+ 




+ 




+ 


1 


+ 


[2,2] 


+ 




+ 




+ 


1 


1 



Table 12: Numbering functions with local properties for / G £i i, / G £1 2, / G >Ci 3, 



g.o.p. 


N=5 


N=5 


N=6 


N=6 


N=7 


N=7 


Total number 




259 




950 




387 


[1] 


95 


+ 


340 


+ 


1,193 


+ 


[1,1] 


50 


+ 


174 


+ 


600 


+ 


[1, 1, 1] 


16 


+ 


58 


+ 


204 


+ 


[1,1,1,1] 


4 


+ 


16 


+ 


60 


+ 


[1,1,1,1,1] 


1 


+ 


4 


+ 


16 


+ 


[1. 1. 1. 1. 1.1] 




+ 


1 


+ 


4 


+ 


[1,1,1,1,1,1,1] 




+ 




+ 


1 


+ 


[2] 


70 


70 


264 


264 


952 


952 


[2,1] 


12 


18 


45 


70 


166 


264 


[1,2] 


6 


+ 


25 


+ 


98 


+ 


[2,2] 


4 


4 


18 


18 


70 


70 


[2,2,1] 


1 


1 


4 


4 


17 


18 


[1,2,2] 




+ 




+ 


1 


+ 


[2,1,2] 




+ 




+ 




+ 


[2,2,2] 




+ 


1 


1 


4 


4 


[2,2,2,1] 




+ 




+ 


1 


1 



Table 13: Numbering functions with local properties for / G £1,5, / G £1,6, / G £1,7. 
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g.o.p. 


N=8 


N=8 


N=9 


N=9 


N=10 


N=10 


Total number 




11,814 




40,503 




13,6946 


[1] 


4,116 


+ 


14,001 


+ 


47,064 


+ 


[1,1] 


2,038 


+ 


6,852 


+ 


22,806 


+ 


[1,1,1] 


700 


+ 


2,366 


+ 


7,896 


+ 


[1,1,1,1] 


214 


+ 


742 


+ 


2,520 


+ 


[1,1,1,1,1] 


60 


+ 


216 


+ 


754 


+ 


[1,1,1,1,1,1] 


16 


+ 


60 


+ 


216 


+ 


[1,1,1,1,1,1,1] 


4 


+ 


16 


+ 


60 


+ 


[1,1,1,1,1,1,1,1] 


1 


+ 


4 


+ 


16 


+ 


[1,1,1,1,1,1,1,1,1] 


- 


+ 


1 


+ 


4 


+ 


[1,1,1,1,1,1,1,1,1,1] 


- 


+ 


- 


+ 


1 


+ 


[2] 


3,356 


3,356 


11,580 


11,580 


39,364 


39,364 


[2,1] 


590 


952 


2,062 


3,356 


7,072 


11,580 


[1,2] 


362 


+ 


1,294 


+ 


4,508 


+ 


[2,2] 


264 


264 


952 


952 


3,356 


3,356 


[2,2,1] 


62 


70 


222 


264 


770 


952 


[1,2,2] 


6 


+ 


28 


+ 


113 


+ 


[2,1,2] 


2 


+ 


14 


+ 


69 


+ 


[2,2,2] 


18 


18 


70 


70 


264 


264 


[2,2,2,1] 


4 


4 


18 


18 


69 


70 


[1,2,2,2] 




+ 




+ 


1 


+ 


[2,1,2,2] 




+ 




+ 




+ 


[2,2,1,2] 




+ 




+ 




+ 


[2,2,2,2] 


1 


1 


4 


4 


18 


18 


[2,2,2,2,1] 




+ 


1 


1 


4 


4 


[1,2,2,2,2] 




+ 




+ 




+ 


[2,1,2,2,2] 




+ 




+ 




+ 


[2,2,1,2,2] 




+ 




+ 




+ 


[2,2,2,1,2] 




+ 




+ 




+ 


[2,2,2,2,2] 




+ 




+ 


1 


1 



Table 14: Numbering functions with local properties for / e jCi^g, f e jCi^g, f e >Ci_io. 
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g.o.p. 


N=ll 


N=ll 


N=12 


N=12 


N=13 


N=13 


Total number 




457,795 




1,515,926 




4,979,777 


[1] 


156,629 


+ 


516,844 


+ 


1,693,073 


+ 


[1,1] 


75,292 


+ 


246,762 


+ 


803,706 


+ 


[1,1,1] 


26,098 


+ 


85,556 


+ 


278,580 


+ 


[1,1.1,1] 


8,434 


+ 


27,904 


+ 


91,488 


+ 


[1,1,1,1,1] 


2,756 


+ 


8,658 


+ 


28,738 


+ 


[1,1,1,1,1,1] 


756 


+ 


2,590 


+ 


8,730 


+ 


[1,1,1,1,1,1,1] 


216 


+ 


756 


+ 


2,592 


+ 


[1,1,1,1,1,1,1,1] 


60 


+ 


216 


+ 


756 


+ 


[1,1,1,1,1,1,1,1,1] 


16 


+ 


60 


+ 


216 


+ 


[1,1,1,1,1,1,1,1,1,1] 


4 


+ 


16 


+ 


60 


+ 


[1,1,1,1,1,1,1,1,1,1,1] 


1 


+ 


4 


+ 


16 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1] 


- 


+ 


1 


+ 


4 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1,1] 


- 


+ 


- 


+ 


1 


+ 


[2] 


132,104 


132,104 


438,846 


438,846 


1,445,258 


1,445,258 


[2,1] 


23,941 


39,364 


80,108 


132,104 


265,548 


438,846 


[1,2] 


15,423 


+ 


51,996 


+ 


173,298 


+ 


[2,2] 


11,580 


11,580 


39,364 


39,364 


132,104 


132,104 


[2,2,1] 


2,634 


3,356 


8,883 


11,580 


29,659 


39,364 


[1,2,2] 


429 


+ 


1,555 


+ 


5,478 


+ 


[2,1,2] 


293 


+ 


1,142 


+ 


4,227 


+ 


[2,2,2] 


952 


952 


3,356 


3,356 


11,580 


11,580 


[2,2,2,1] 


255 


264 


899 


952 


3,098 


3,356 


[1,2,2,2] 


7 


+ 


35 


+ 


152 


+ 


[2,1,2,2] 


2 


+ 


16 


+ 


86 


+ 


[2,2,1,2] 


- 


+ 


2 


+ 


20 


+ 


[2,2,2,2] 


70 


70 


264 


264 


952 


952 


[2,2,2,2,1] 


18 


18 


70 


70 


263 


264 


[1,2,2,2,2] 


- 


+ 


- 


+ 


1 


+ 


[2,1,2,2,2] 


- 


+ 


- 


+ 


- 


+ 


[2,2,1,2,2] 




+ 




+ 




+ 


[2,2,2,1,2] 




+ 




+ 




+ 


[2, 2, 2, 2, 2] 


4 


4 


18 


18 


70 


70 


[2,2,2,2,2,1] 


1 


1 


4 


4 


18 


18 


[1,2,2,2,2,2] 




+ 




+ 




+ 


[2,2,2,2,2,2] 




+ 


1 


1 


4 


4 


[2,2,2,2,2,2,1] 




+ 




+ 


1 


1 



Table 15: Numbering functions with local properties for / e jCm, f e jCiu, 
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g.o.p. 


N=14 


N=14 


N=15 


N=15 


Total number 




16,246,924 




52,694,573 


[1] 


5,511,218 


+ 


17,841,247 


+ 


[1,1] 


2,603,258 


+ 


8,391,360 


+ 


[1,1,1] 


901,802 


+ 


2,904,592 


+ 


[1,1,1,1] 


297,728 


+ 


962,888 


+ 


[1,1,1,1,1] 


94,440 


+ 


307,848 


+ 


[1,1,1,1,1,1] 


29,050 


+ 


95,676 


+ 


[1,1,1,1,1,1,1] 


8,746 


+ 


29,140 


+ 


[1,1,1,1,1,1,1,1] 


2,592 


+ 


8,748 


+ 


[1,1,1,1,1,1,1,1,1] 


756 


+ 


2,592 


+ 


[1,1,1,1,1,1,1,1,1,1] 


216 


+ 


756 


+ 


[1,1,1,1,1,1,1,1,1,1,1] 


60 


+ 


216 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1] 


16 


+ 


60 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1,1] 


4 


+ 


16 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1,1,1] 


1 


+ 


4 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1] 


- 


+ 


1 


+ 


[2] 


4,725,220 


4,725,220 


15,352,392 


15,352,392 


[2,1] 


873,149 


1,445,258 


2,851,350 


+ 


[1,2] 


572,109 


+ 


1,873,870 


+ 


[2,2] 


438,846 


438,846 


1,445,258 


1,445,258 


[2,2,1] 


98,135 


132,104 


322,310 


438,846 


[1,2,2] 


18,873 


+ 


63,967 


+ 


[2,1,2] 


15,096 


+ 


52,569 


+ 


[2,2,2] 


39,364 


39,364 


132,104 


132,104 


[2,2,2,1] 


10,460 


11,580 


34,845 


39,364 


[1,2,2,2] 


605 


+ 


2,282 


+ 


[2,1,2,2] 


389 


+ 


1,596 


+ 


[2, 2, 1, 2] 


126 


+ 


641 


+ 


[2,2,2,2] 


3,356 


3,356 


11,580 


11,580 


[2,2,2,2,1] 


942 


952 


3,292 


3,356 


[1,2,2,2,2] 


8 


+ 


44 


+ 


[2,1,2,2,2] 


2 


+ 


18 


+ 


[2,2,1,2,2] 


- 


+ 


2 


+ 


[2,2,2,1,2] 


- 


+ 


- 


+ 


[2,2,2,2,2] 


264 


264 


952 


952 


[2,2,2,2,2,1] 


70 


70 


264 


264 


[1,2,2,2,2,2] 




+ 




+ 


[2,2,2,2,2,2] 


18 


18 


70 


70 


[2,2,2,2,2,2,1] 


4 


4 


18 


18 


[2,2,2,2,2,2,2] 


1 


1 


4 


4 


[2,2,2,2,2,2,2,1] 




+ 


1 


1 



Table 16: Numbering functions with local properties for / e jCi^u, f G £1,15. 
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g.o.p. 


N=16 


N=16 


Total number 




170,028,792 


[1] 


57,477,542 


+ 


[1,1] 


26,932,398 


+ 


[1,1,1] 


9,314,088 


+ 


[1,1,1,1] 


3,097,650 


+ 


[1,1,1,1,1] 


996,764 


+ 


[1,1,1,1,1,1] 


312,456 


+ 


[1,1,1,1,1,1,1] 


96,096 


+ 


[1,1,1,1,1,1,1,1] 


29,158 


+ 


[1,1,1,1,1,1,1,1,1] 


8,748 


+ 


[1,1,1,1,1,1,1,1,1,1] 


2,592 


+ 


[1,1,1,1,1,1,1,1,1,1,1] 


756 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1] 


216 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1,1] 


60 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1,1,1] 


16 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1] 


4 


+ 


[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1] 


1 


+ 


[2] 


49,610,818 


49,610,818 


[2,1] 


9,255,822 


15,352,392 


[1,2] 


6,096,570 


+ 


[2,2] 


4,725,220 


4,725,220 


[2,2,1] 


1,051,686 


1,445,258 


[1,2,2] 


213,975 


+ 


[2,1,2] 


179,597 


+ 


[2,2,2] 


438,846 


438,846 


[2,2,2,1] 


114,798 


132,104 


[1,2,2,2] 


8,284 


+ 


[2,1,2,2] 


6,146 


+ 


[2,2,1,2] 


2,876 


+ 


[2,2,2,2] 


39,364 


39,364 


[2,2,2,2,1] 


11,246 


11,580 


[1,2,2,2,2] 


204 


+ 


[2,1,2,2,2] 


106 


+ 


[2,2,1,2,2] 


22 


+ 


[2,2,2,1,2] 


2 


+ 


[2,2,2,2,2] 


3,356 


3,356 


[2,2,2,2,2,1] 


951 


952 


[1,2,2,2,2,2] 


1 


+ 


[2,2,2,2,2,2] 


264 


264 


[2,2,2,2,2,2, 1] 


70 


70 


[2,2,2,2,2,2,2] 


18 


18 


[2,2,2,2,2,2,2, 1] 


4 


4 


[2,2,2,2,2,2,2,2] 


1 


1 



Table 17: Numbering functions with local properties for / e Ci^ 
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maximal period 


modulus 


gop number 


functions number 
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9,992 
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41 
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Q 
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51 


o 
O 
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y 


1 00 01 
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OO 


4 


5 


19 
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DU 
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b 


zo 
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4 


ID 


o7 
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oz 


4 




44 
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oi 


D 
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OZ 


D 
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iU4 
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CO 
OO 


( 
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QQ 
OO 
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1 n 
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10 AQK KCiP. 

lo,4o0,0Ub 
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oy 


Q 

o 


1 n 




14,Dy2,D0o 


on 


Q 

o 


1 n 




lo,yo4, /o2 


no 

yz 


o 
o 


11) 


1 OO 

ioz 


lo, ( /o,2o4 


yo 


Q 

o 




1 QQ 


on 0^0 no /I 
2U,202,Uo4 




Q 

o 


1 

-LU 






95 


8 


10 


195 


23,021,112 


96 


8 


10 


242 


24,479,312 


97 


8 


10 


298 


26,163,582 


98 


8 


10 


335 


28,285,274 


99 


9 


10 


377 


30,861,396 


100 


10 


10 


463 


34,086,310 


101 


10 


10 


484 


37,553,504 



18: Numerical study of the set Cat.q.N for N — 10, t — 5, ai — 20, a2 
5, q;4 = 3 and = 1, for q = 35, 101 
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5.2 Orbits and gop of locally bounded range function sets 

Consider now the set : 

^ctt,,,N = {f such that Vp, < p < iV - r - 1, E - /(p + r)| < 

r=t 

q} f]{f G J-'n such that Wp, t < p < N — (^r\f{p) — f{p — r)\ < q} for the vector 

r=l 

Ot = (ai, 0(2, ■ ■ ■ , ctt) G N*, for q eN. 

The functions belonging to these sets show a kind of "rigidity": the less is q, the 
more "rigid" is the function. Furthermore, the maximal length of a periodic orbit 
increases with q, and so the number of gop Wi'^c^gN) the maximal modulus 
of the gop. As an example, we explore numerically the case: iV = 10, t = 5, 
ai = 20, 0^2 = 10, = 5, 0^4 = 3 and = 1, for g = 35, . . . , 101. The results are 
displayed in table [181 In this table "modulus" means the maximal modulus of the 
gop belonging to this set for the corresponding value of q in the row, "gop number" 
stands for Wi'^al g j^) "functions number" for tl-^^q^- One can point out 
that for the particular function at of the example; it is possible to find 10 intervals 
Ji, /2, . . . , /lo C N such that if g G /r then there is no periodic orbit whose period is 
strictly greater than r, (e.g., Iq = |70,82]). Furthermore it is possible to split these 
intervals into subintervals Ir^s in which ^Q{C-^ ^ ^) is constant when q thumbs Ir^s- 
This is not the case for tj-^;^ 



6 Conclusion 



The behaviour of a discrete dynamical system associated to a function on finite 
ordered set X is not easily predictable. Such a system can only exhibit periodic 
orbits. All the orbits can be explicitly computed; all together they form a global 
orbit pattern. We formalise such a gop as the ordered set of periods when the 
initial value thumbs X in the increasing order. We are able to predict, using closed 
formulas, the number of gop for the set J-'n of all the functions on X. We introduce 
special subsets of J-'jy in order to understand more precisely the behaviour of the 
dynamical system. We explore by computer experiments theses sets, they show 
interesting patterns of gop. Further study is needed to understand the behaviour of 
dynamical systems associated to functions belonging to these sets. 
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